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p ■ Abstract 

OO . The heavy jet mass distribution in e~^e~ colhsions is computed to next-to-next-to- 

next-to leading logarithmic (N^LL) and next-to-next-to leading fixed order accuracy 

i-Q ' (NNLO). The singular terms predicted from the resummed distribution are confirmed 

by the fixed order distributions allowing a precise extraction of the unknown soft function 

^h! coefficients. A number of quantitative and qualitative comparisons of heavy jet mass and 

r"! I the related thrust distribution are made. From fitting to ALEPH data, a value of a^ 

is extracted, as{mz) = 0.1220 it 0.0031, which is larger than, but not in confiict with, 

^T) . the corresponding value for thrust. A weighted average of the two produces as{mz) = 

K*" I 0.1193 zb 0.0027, consistent with the world average. A study of the non-perturbative 

^L ' corrections shows that the flat direction observed for thrust between a^ and a simple 

\^ . non-perturbative shape parameter is not lifted in combining with heavy jet mass. The 



Monte Carlo treatment of hadronization gives qualitatively different results for thrust 
VQ ' and heavy jet mass, and we conclude that it cannot be trusted to add power corrections 

to the event shape distributions at this accuracy. Whether a more sophisticated effective 
field theory approach to power corrections can reconcile the thrust and heavy jet mass 
distributions remains an open question. 
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1 Introduction 

Event shapes in e~^e~ collisions provide some of the best ways to test QCD and the standard 
model. At high energies, where QCD is perturbative, event shapes lead to some of the world's 
most precise measurements of the strong coupling constant a^. Recently, a number of theoret- 
ical advances have led to renewed interest in event shapes and the a^ measurements. First, the 
NNLO fixed order Feynman diagrams were calculated [H [21 El S] • This allowed the prediction 
of all event shapes to order a^. Second, advances [5l [6l [7] in Soft-CoUinear Effective Theory 
(SCET) in m [To] have allowed resummation of the large logarithmic corrections to thrust to 
N'^LL accuracy [11]. Previous calculations were at NLO [12] and NLL [13]. Very recently, 
a full effective field theory analysis of a single event shape, thrust, has been completed [H], 
including additionally non-perturbative considerations. The resulting a^ extraction is com- 
petitive with the PDG world average [15], however it differs significantly from measurements 
using lattice QCD and r decays (see [16] for a review). 

Although the as measurement with thrust is extremely precise, there are many reasons to 
study additional event shapes as well. The main advantage is that there may be systematic 
effects in a single event shape pulling Og in a certain direction, which are not universal. In 
fact, as observed in [17] from an NLL-I-NNLO analysis, there seem to be two classes of event 
shapes, the first including thrust, the C-parameter and total jet broadening, while the second 
includes heavy jet mass, wide jet broadening and the two-to-three jet transition parameter 
1/3. The values of a^ extracted from the two classes at NLL+NNLO tend to have around a 
5% systematic difference which the authors attribute to missing higher order corrections. In a 
recent world average of a^ [TB], the thrust measurement using SCET was not included because 
of a concern over precisely this kind of systematic uncertainty. In this paper, we correct that 
concern with a N^LL calculation of an event shape from the second class, heavy jet mass. 

In addition to being useful for measuring as, the heavy jet mass distribution allows us to 
explore other aspects of resummation. Indeed, there are few hadronic observables which have 
been calculated this accurately and for which there is data. Heavy jet mass involves a soft 
function which cannot be written in terms of only a single scale. These types of soft functions 
promise to play an important role in resummation at hadron colliders and only beginning 
to be explored [191 [201 121] . We will discuss constraints on the soft function, and perform a 
numerical study of the parts that are not known, similar to what was done in [TT] and [22] • 

Hadronization is another issue which having a second event shape may help understand. 
In the fit to a^ with the thrust distribution [11] , it was observed that a decrease in as could be 
compensated for with a single non-perturbative parameter with only a small effect on the x^ 
of the fit. Having another observable for which the same non-perturbative parameters can be 
fit can possibly remove this flat direction. The hadronization issue is also important for Monte 
Carlo simulations. With a more accurate theoretical calculations, we can explore whether the 
approximations in PYTHIA [23] allow for an adequate description of thrust and heavy jet mass 
simultaneously. 

As a brief outline of our findings, we begin in Section [2] with an overview of the SCET 
distributions. The hemisphere soft function is studied and singular terms in the heavy jet 
mass distribution are compared to the fixed order calculation in Section [31 We found a 
mild inconsistency with the analytic results from SCET and the numerical calculations of the 



NLO and NNLO distributions. After completing the original study, we were able to resolve 
this inconsistency, which was due to Monte Carlo convergence problems, by taking a very 
low numerical infrared cutoff, as discussed in a note added at the end of this paper. In 
Section H] we fit for a^. The fit for heavy jet mass to the lep data from aleph [23] leads to 
as{mz) = 0.1220 ± 0.0031. This value is higher than the value from thrust using exactly the 
same technique, as{Taz) = 0.1175 ±0.0026. Assuming 100% correlation gives an average value 
of asirriz) = 0.1193 ± 0.0027 which is very close to the recent average in [TB]. We also find 
that convergence of the perturbation series for heavy jet mass with resummation is, like for 
thrust, significantly better than the convergence of the fixed order calculation. In Section |5l a 
comparison of the data to pythia shows that while pythia agrees with the thrust data almost 
perfectly, it has trouble matching the heavy jet mass distribution. Moreover, the hadronization 
corrections in PYTHIA move the curve in the wrong direction for heavy jet mass. Concluding 
that the Monte Carlo hadronization model is incompatible with the high precision theoretical 
calculation, we explore non-perturbative corrections in SCET with a simple shape function. 
We find that to the order we are working, this simple shape function cannot simultaneously 
describe the thrust and heavy jet mass distributions. We provide an expanded summary, 
discussion and comparison to previous results in Section O 

2 Thrust and Heavy Jet Mass in SCET 

Thrust and heavy jet mass are defined as follows. One first finds the thrust axis, through 

T = max ' — — , (1) 

l^ilPil 

where the sum is over all momentum 3- vectors pj in the event, and the maximum is over all 
unit 3- vectors n. We use r = 1 — T to measure thrust. Once the thrust axis is known, it can 
be used to split the event into two hemispheres. We define P^ and Pj^ to be the four momenta 
of the sum of all the radiation going into each hemisphere, and Ml = \/Pl and Mji = \/P^ 
to be the hemisphere masses. Heavy jet mass p is defined as the larger of the two hemisphere 
masses squared, normalized to the the center of mass energy Q, 

p^i_inax(Ml,M^). (2) 

When p is small, r is also small, both hemisphere masses are small, and the event appears 
to have two back-to-back pencil-like jets. In this threshold limit, the thrust axis aligns with 
the jet axis and Q^t approaches the sum of the two hemisphere masses squared M^ + M^ = 

gv + o(r2). 

It follows that both thrust, up to power corrections, and heavy jet mass can be written as 
integrals over the doubly differential hemisphere mass distribution: 



In SCET, the doubly differential hemisphere mass distribution is calculable in the threshold 
limit. The appropriate factorization theorem in SCET was first derived in [5] for the related 
process of ti production. This theorem was then shown to allow for the calculation of event 
shapes in |6], where matched and resummed thrust and heavy jet mass distributions in SCET 
were first presented. Previously, resummation of heavy jet mass was only possible at NLL 
accuracy [13]. The first event shape resummed to N^LL was thrust, in [11]. Monte Carlo 
based hadronization corrections were included in [25] to produce a strong model-independent 
gluino mass bound. Recently, power corrections for thrust were studied within the effective 
field theory approach in [2]. 

The factorization theorem allows us to write the hemisphere mass distribution as 
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H{Q\ fi) jdkLdkR J{Ml - QkL, /i) J{M'^ - Qkn, fi) S{kL, kn, fi) . (5) 



The subscript on a2 is a reminder that this expression holds in the two-jet region. Here, 
H[Q'^,fi) is the hard function. It is calculated in matching SCET to QCD and contains 
information about the modes of QCD that are not in SCET. J{p'^,fi) is a jet function. It is 
derived in a matching calculation from a theory with soft and coUinear modes to a theory 
with just soft modes. The hard function was calculated in [26] and the jet functions in [7]. 
Finally, S^k^, kptyH) is the hemisphere soft function which is derived from integrating out the 
remaining soft modes. 

The doubly differential hemisphere mass distribution, Eq. ([5]), is observable, and therefore 
must be independent of renormalization group scale fi. Demanding /^-independence leads 
to a renormalization group equation which is easiest to express in Laplace space, where the 
convolutions turn into products. The Laplace transform is defined by 

f{uL, ur) = I dMldMJ,e-^-''le-^-'''-f{Ml Ml) (6) 

which can be applied to the differential cross section and to the jet and soft functions separately. 
We generally express the Laplace transformed distributions as functions of Li = hi{^VLe""^) 
and L2 = ln(/iz/^e'''^). Since the entire /i-dependence of the hard and jet functions is known, 
the /i-dependence of the soft function is completely fixed by renormalization group invariance 
(see [TT] for more details). The result is that the hemisphere soft function itself factorizes into 
the form [3 El [23 |22] 

s(Li,L2,/i) = s^(Li,/i)s^(L2,/i)s/(Li - L2) , (7) 

where all the /i-dependence is contained in the function 's^{L, fi) which is known to N'^LL 
accuracy. Since Li — L2 = ln(i^i/z//j), the function 'sf{Li — L2) is /i-independent. We discuss 
the soft function more in the next section. 

Putting together the hard and jet functions with the soft function written in this way 
produces an analytic expression for the doubly differential jet mass distribution. For thrust, 
the result is [TT] 
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with ?7 = 4:Ar{fij, fig) and the thrust soft function st{L,^) is defined by 



(9) 



The definitions of the RG kernels A^iiy, /i) and ^(z/, /i) as well as the fixed order hard and jet 
functions, H{Q'^,jj) and j{L,fi) and their anomalous dimensions can be found in [H]. Note 
that only one value of the unknown function S/(L) is required for thrust. 
For heavy jet mass, the distribution is similar 
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where ?7i = ?72 = 2Ar(/ij, ps). In contrast to thrust, for heavy jet mass the full functional form 
of 'sf{L) is needed. For N^LL precision, we need to know the hemisphere soft function, and 
hence 'sf{L) to two-loop order (NLO). Actually, to this order, we only need one projection 
of the hemisphere soft function. For three-loop matching (NNLO), we need an additional 
projection. These projections will be discussed in the next section. 

One interesting feature of the hemisphere mass distribution is that the soft interference 
effects in 'sf{L) are only relevant at a^, which is appropriate for N^LL resummation. Up 
to NNLL accuracy, the doubly differential distribution is simply the product of the mass 
distributions in the two hemispheres. Explicitly, 
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and 1] = 2Ar{fij, fig)- Since, for NNLL resummation, the hard and jet functions are only 
needed to 0{as), the square-roots above simply mean take one half of the «s pieces. The fact 
that the distribution splits up in this way was observed at NLL level in [6], and is essential to 
the traditional NLL resummation [13]. This simplified factorization suggests that it may be 
possible to calculate observables involving many more jets with NNLL resummation without 
having to disentangle soft interference effects. Note that this factorization does not guarantee 
that large logs of Ml/Mj^ can be resummed. However, it is possible that the calculation of 
observables with only one scale, such as the sum of many jet masses, or a maximal jet mass, 
will simplify with SCET. 

3 Hemisphere Soft Function and Comparison to Fixed 
Order 

The hemisphere soft function has been studied briefly in [5l[6l[27] and more thoroughly in [22] . 
It is a function of two scales, k^ and fc/j as well as the renormalization group scale /i. If n^ 
is the direction of the left hemisphere and k^ is the sum of the momenta of all the soft 
radiation entering this hemisphere, then ki is the component of /c^ backwards to n^. That 
is ki = {kL-riL). kji is defined analogously. The soft function can be factorized into a 
perturbative, partonic part, and non-perturbative contribution which has support of order 
Aqcd- For now we deal only with the perturbative part, discussing non-perturbative effects 
in Section 

As we have noted, the soft function itself factorizes. 

s(Li,L2,/i) = s^(Li,/i)s^(L2,/i)s/(Li -La) (13) 

where 's{Li, -^2, /i) is the Laplace transform of S^ki, kji, /^), as in Eq. ([6]), and Li = ln(/iz/ie'^^), 
L2 = ln(/iz/ije'''^). The function s^(L,/i) is completely fixed by RG invariance in terms of the 
hard and jet anomalous dimensions. It can be calculated in perturbation theory by demanding 
Eq. (fTOj) be independent of fj,. This gives 



s,(L,/i) = exp [ (^) {-L'To + L7o^) + (^)' Ql^^/^oFo + L\-T^ - /3o7o^) + ^(7f )) 
+ {^y (-^^'/3oro + Il'WiTo + 2/3ori + 2/3o%o^) + L'{-T2 - ^il - 2/3o7f ) + L{^l)\ 

■ (14) 

The yU-independent part S'/(L) must satisfy a number of constraints, as discussed in [22] . 

First of all, since the soft function is symmetric in the two hemispheres, ?/(L) must be an 
even function of L. Second of all, we know the function to order as by explicit calculation. 



Writing 

the one-loop result is that 

Sf,{L) = -CFn\ (16) 

The authors of [22] also observed that 'sf{L) is constrained by the non-Abelian exponentiation 
theorem. Non-Abelian exponentiation implies constraints on powers of logarithms of fi in 
the full soft function. These constraints are satisfied by the explicit solution, since 's^{L,fi) 
is an exponential. The theorem also restricts the Cp color structure in the soft function to 
be completely determined by the one- loop result, Eq. f lTO]) . Beyond this, however, S/(L) is 
unconstrained. It may even have more general dependence on L than logarithms. To determine 
'sf{L), we must calculate the soft function perturbatively. The one- loop calculation has been 
done but the two-loop calculation, which is required for N'^LL resummation, has not. 

A simple alternative to calculating 'sf{L) at NNLO is to extract projections of s'/(L) from 
numerical comparisons to event shape calculations in full QCD. For example, thrust is only 
sensitive to sj(0). Writing 
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and comparing to Eq. ( fT6l) . we see that cf = —Cp'n'^. The two-loop constant was determined 
numerically in [H] with the use of the EVENT 2 program [28] . The result is 

cf = (58 ± 2)C| + (-60 ± l)CpC^ + (43 ± l)Cpr^n/ (Becher and Schwartz) (18) 



This is in conflict with the prediction from non-Abelian exponentiation, which requires the 
C'l' factor be ^n^Cp = 48.70 p. The two-loop constant was also determined in [22], using the 
same technique but imposing non-Abelian exponentiation. They found 
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Cj, + (-59 ± 2)CfCa + (44 ± 3)CFTFnf (Hoang and Kluth) 



(19) 



The two results agree, except for the Cp term. Indeed, the Cp term seems to indicate a 
disagreement between the numerical results of the event 2 program and the prediction from 
non-Abelian exponentiation. Since the uncertainty in Eq ( ITSl) is too small to explain this 
disagreement, it is reasonable also to expect the other color structures to be off. We should 
therefore allow for a systematic uncertainty on these flts in addition to what is presented, 
which is essentially a statistical uncertainty associated with the flt. We discuss this more 
below. 

Event shapes other than thrust are sensitive to the form of s'/(L), not just s'j(O). This can 
be seen, for example, by the form of the heavy jet mass distribution in Eq. (TTUj) . For N'^LL 
resummation, the flxed order expansion is required to a^. The contribution at this order 
involving 'sf{L) requires at most 'sf2{L), with the jet and hard functions at their tree-level 
values. Thus, the required projection of the s'j(L) for heavy jet mass is 
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Figure 1: A comparison of the full fixed-order calculations and expanded SCET at NLO. 
Shown is the difference -^pA^ = p(B{p) — Db{p)), where B{p) is the full NLO 5-function, 
calculated with event 2 and Db{p) is the singular part, calculated with SCET. The differences 
are separated by color structure, with the sum also shown. The kink at p = | is the maximum 
heavy jet mass for a 3-particle final state. (See also Figure [TTl) 



The integral representation of cf is suggestive of a deeper relation between heavy jet mass 
and the hemisphere mass distribution, however we do not have a physical explanation of why 
this particular moment appears. If S/(L) is a polynomial, this moment is very simple. For 
example, if we assume 



ML) = 1 + (^) cf + (^)' [cf + 4,L= + c,|,L'] 



then 
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At NLO, the singular part of the heavy jet mass distribution only depends on S'/(L) through 
cf . Thus, we can fit cf numerically the same way cf is fit with thrust. 

To determine cf we use the same technique used in [11] for thrust, and in [22] for a one- 



-2p 

parameter family of event shapes. The basic idea is that the singular part of the heavy jet 
mass distribution is known analytically, through SCET. The difference between the exact NLO 
heavy jet mass distribution and this singular part is finite and can be integrated numerically. 
This integral is then the total cross section at NLO minus the integral of the singular part 
which is calculable analytically and depends on the constants cf for thrust and cf for heavy 
jet mass. 
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Figure 2: Extraction of the two-loop constants in the soft function. The points correspond 
to the value of a lower bound pmin applied to the fixed-order calculation. The lines are 
interpolations among the points from pmin = 0.002 to pmin = 0.005 extrapolated to p = to 
extract the constants. (See also Figure [T2l) 



In more detail, the exact heavy jet mass distribution can be expanded as a series in as 
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Each term in this series is singular at p 
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0. The singular parts can be written as a sum of 
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The functions Ds, Da{p), Db{p), and Dc{p) are calculable in SCET and we give them in 
Appendix O Up to order a^, the only dependence on the unknown soft function coefficient 
cf is in Ds, thus the shape of the singular part of the NLO distribution is known completely. 
The corresponding exact distributions in perturbative QCD have been calculated for p > 
analytically for the A function, and numerically for the B and C functions. Since SCET 
produces the entire singular part of the distributions, the combination 
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should vanish at p = 0. We show this difference separated by color structure in Figure [TJ The 
B functions are calculated using the Monte Carlo program event 2 [28] with 10^° events. 
Curiously, while the Cp and C^ color structures do not seem to go to as p — )■ 0, their sum 
does. 
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Figure 3: Comparison of the full NNLO heavy jet mass distribution C{p) (dashed black 
histograms) [22] and the singular terms Dc{p) (blue curves). The light-red areas are an 
estimate of the statistical uncertainty. The uncertainties on cf and c|^, in Eq fl33|) . are not 
visible. The disagreement at very small p is due to the infrared cut-off of hq = 10~'^ for 
the NNLO calculation. It is expected that the agreement would improve if this cutoff were 
lowered, as can be seen in the analogous thrust plot in [11]. (See also Figure | 



With these difference functions, it is straightforward to extract cf^ as in [TT] and [22] . 
Although the difference pB{p) — pDsip) is regular at p = 0, the two functions are separately 
divergent. Since Db{p) is only known numerically, the difference is numerically unstable at 
small p. To do the cf extraction, we use the same procedure as in [11] and impose an lower 
bound and take the limit that the bound is removed. The extracted values as a function of this 
lower bound pmin are shown in Figure |2J We then fit a line in the region 0.002 < pmin < 0.005 
and extrapolate to pmin = 0. The result is 



cf = (58 ± 2)C| + (-41 ± 2)CfCa + (50 ± l)CpT^n/ . 



(26) 



Note that the Cp and Ca curves have problems at small Pmim in agreement with what is seen 
in Figure [H Since the approach is linear up to around p^i^ ~ 0.002, it is likely that this 
divergence is an unphysical systematic problem with the Monte Carlo, and not due to low 
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Figure 4: The difference, ^pA^ = p{C{p) — Dc{p)), between the full NNLO heavy jet mass 
distribution and the singular terms. The light-red areas are an estimate of the statistical 
uncertainty, from [29]. The blue band is the uncertainty due to cf and c|^, in Eq ( 133|1 . These 
curves should all go to at p = 0. The C{p) distributions are all calculated with an infrared 
cutoff of Ho = 10"^. (See Figure [H] for the same figure with i/q = 10"'''.) 

statistics or a discrepancy with theory. \j As with cf , non-Abelian exponentiation implies that 
the Cp term should be ^n^ 



48.7. Thus, our uncertainty of cf from the extrapolation is 



probably too small and we will therefore inflate the errors by a factor of 5. (See also Eq. ( 1501) 
and Figures ITT] and fT2] for an update.) 

To calculate the heavy jet mass distribution to N'^LL-|-NNLO accuracy, we must match to 
the NNLO fixed order distribution. This requires the singular parts of heavy jet mass to a^, 
that is, the function Dc{p) in Eq. flMl) . To derive this, we do not need the finite part of the 
soft function at a^, 'sf^{L), since this piece only contributes to the a^ part of Ds-, which is not 
required for matching. We do, however, need another projection of the a^ soft function, of 
the form 



cfc = i.dvi + 9n2)sf2idn, - d, 
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^We thank A. Hoang for a discussion of this point. 
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Again, we have no physical explanation of the intriguing integral definition in the second line. 
This projection also simplifies with a polynomial soft function. For example, with Eq (l2Tll 



c: 
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4C3cfi + (-Stt^Cs + 48C5)cfQ . (29) 

The prediction from SCET for Dc{p) with its explicit dependence only on cf and cf^ is given 
in Appendix O There are only three color structures which depend on cf^ at all. 

In order to extract the L dependence of s/2(i^), we could attempt to fit cf^ with the shapes 
of the NNLO distributions. An alternative, as pursued by Hoang and Kluth in [22], is to use 
the other event shapes beyond thrust and heavy jet mass at NLO. These authors considered 
a weighted sum of the jet masses, Tq = j^{aMl + Mj^/Q"^. This form leads to a singular 
distribution which depends on sj2(lna), hence combining event shapes with different a can 
probe the entire function 'sf2{L). Their fits show good agreement with the form 

'^fW-l+{^)cU{^y[cU4LL']. (30) 

which they have argued is likely to be the exact 2-loop soft function. We will therefore assume 
this form of the soft function as well, in order to proceed with the N^LL+NNLO a^ fits, g 

With this soft function and the thrust fit values in Eq. (ITH]) . our fit for cf^ translates into 
a fit for cfi (cf. Eq.(l22]) with cfg = 0) The result is 

4l = (0 ± 2)C| + (-5.8 ± 1.5)CfCa + (-2.2 ± l)Ci.Ti.n/ . (31) 

Using a similar technique, but imposing the constraint from non-Abelian exponentiation, 
Hoang and Kluth found results consistent with ours 



c 



f^ = (O)C^ + (-6.5 ± 2)CfCa + (1.3 ± 2)CFTFnf (Hoang and Kluth) (32) 



Note that for cf^, the Cp coefficient comes out to be consistent with the prediction from non- 
Abelian exponentiation. Since cf^ comes from the difference between the values extracted from 
thrust and the values extracted from heavy jet mass, the systematic problem with event 2 
may be cancelling in the difference. Thus, we will inflate our uncertainties on cfj, by only a 
factor of 2. 

In summary, for the rest of this paper, we will take 

ML)-l^{l^)c^,+{^y[4 + 4LL'-] (33) 

cf = —Cl + (-60 ± 10)CfCa + (43 ± 5)CFTFnf (34) 



^ There is a subtlety about these Tq event shapes because of non-global logarithms [30]. For example, for 
very large or small a, these event shapes reduce to the left or right hemisphere mass, which are known to have 
non-global logs. Since r^ — )■ forces the masslcss dijet threshold, in which the SCET factorization theorem 
is derived, only up to corrections of order In a, it is not completely clear that SCET will reproduce all of the 
a-dependence of the singular terms in r^ . 
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4l = (O)C^ + (-6 ± 3)CfCa + (-2 ± 2)CFTFnf , (35) 

so that 

4=4C3cfi and cI = cI-—cIl- (36) 

The uncertainty on a^ due to the uncertainty on these numbers will be included in the fits. 

Before moving on the ag extraction, we can compare the SCET prediction for the singular 
parts of the NNLO distribution to the exact results, as was done for thrust in [TT]. To do this, 
we use Dc{p) from Appendix [0 with the substitutions in Eq. ( l36l) . This lets us compare to the 
C functions in the NNLO distribution, from [22] . Plots of pDc and pC are shown in Figure [3] 
as functions of log p. The uncertainty on cf and cf^ is included, but not visible in these plots. 
Although the agreement is not perfect at very small p, it is expected to improve, as we seen 
for thrust in [11] , as the the infrared cutoff used in the NNLO calculation is reduced from the 
value Ho = 10"^ used here. A version of this plot with cutoff yo = 10"^ has been included as 
Figure [131 confirming our expectations. 

The difference between the full NNLO distribution and its singular parts, as in Eq. (l25l) . is 
shown in Figure IH These curves, for all color structures, should go to zero at p = 0. For most 
of the color structures, this looks plausible, although the 1/A^^ color structure, corresponding 
to the a^Cp coefficient in the heavy jet mass distribution which is fixed by non-Abelian 
exponentiation, looks a bit suspicious. Because this constant is known, we have not included 
an associated uncertainty. The discrepancy is likely due to the infrared cutoff yo = 10~^ used 
for these plots (an update with y^ = 10~^ is included as Figure [14]). Note that even if the 
Ansatz in Eq. (1331) were wrong, a general dependence on cf^ will only affect some of the color 
structures, and even then would only generate an overall up or down shift in these curves (cf. 
the form of Dq in Eq. ([69]) ) . 



4 as extraction and error analysis 

In the previous section, we determined the unknown coefficients in the hemisphere soft function 
and checked the singular terms against the exact NLO and NNLO heavy jet mass distributions. 
Now we are ready to compare to data and fit for the strong coupling constant a^. The 
procedure we follow is identical to the procedure used for thrust in [11] , so we refer the reader 
to that paper for missing details. 

For heavy jet mass, as for thrust, we match to the fixed order distribution via 

r{p) , (37) 







1 dcr 1 d(T2 






(To dp o-Q dp 


with 






r{p) = 


^(i)'-(^)- 


- ^-("'i + (2";) l^*"' - 



' [C(p) - Dcip)] , (38) 

and Da,Db and Dc are given in Appendix [C] The A function is known analytically, and is 
the same as for thrust (see [6]). For B{p) we use the output of event 2 [28], and for C{p) we 
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Figure 5: Relative error for best fit to aleph data at 91.2 GeV. Tlie inner green band includes 
only statistical uncertainty, while the outer yellow band includes statistical, systematic and 
generator uncertainties. The solid black line is for as{Mz) = 0.1214, the best fit value for 
heavy jet mass, while the dashed line has as{'mz) = 0.1168, the best fit for thrust. The fit 
ranges, 0.08 < p < 0.18 and 0.1 < r < 0.24, are taken from [18]. 



use the NNLO calculation which has been provided to us by the authors of [3] . We normalize 



to the total hadronic cross section at order a^, which is 



Chad 



4tt ^ ' V47r 



CfCa 



123 



44C3 + CpTFUf (-22 + I6C3 



Cl- 



■ (39) 



Since the data is binned, what we actually use for the theory prediction is the difference 
between the integrated heavy jet mass distribution evaluated at the bin edges: Rp{p2) — Rp{pi) ■ 
Our fit ranges are chosen to be the same as in [18], so that we can use their values for the 
systematic experimental uncertainties. 

The left panel of Figure [5] shows a comparison of the theory prediction for heavy jet mass 
to the ALEPH data at 91.2 GeV. These curves use the default scale choices 



Idh = Q, Pj = Qy/P, Ps = Qp- 



(40) 



These scales are the natural ones to minimize the large logarithms, and can be read off the 
formula in Eq. (fTOl) . The best fit value of a^ for heavy jet mass is asimz) = 0.1214. We 
show also in the same figure, the heavy jet mass distribution for as{mz) = 0.1168, which is 
the value of a^ derived in [TT] from the fit to the thrust distribution at the same energy. In 
the right panel of Figure [5], we show a comparison to data for thrust, with the same values 
of as- Overall, the fit to thrust is a much better fit. For heavy jet mass, the best fit gives 
X^/d.o.f.=67/9 using statistical uncertainties only, while for thrust, x^/d.o.f.=32.5/13. The 
relatively poor fit for heavy jet mass can be plainly seen in the figure. For thrust, the relative 
distribution is fiat over the fit range (dashed curve, right panel), while for heavy jet mass, 
it is increasing (solid curve, left panel). This coordinates with the relatively larger power 
corrections that we will find in the next section. 

Next, we look at the uncertainties on the theoretical prediction. As with thrust, in jllj . 
we consider first separate variations of ph, pj, ps and the scale pm where the matching is done 
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Figure 6: Convergence of resummed and fixed-order distributions, aleph data (red) at 91.2 
GeV is included for reference. All plots have as{mz) = 0.1214. 

by factors of 2. Figure E] shows the effect of the envelope of these variations on the heavy jet 
mass distribution, for four orders in perturbation theory. We use the same definitions for the 
various orders as in |llj : 



Order 


resum. 


r 

-■- cusp 


7n 


Cn 


matching 


P* order 


NIL 


2-loop 


1-loop 


tree 


- 


2'^^^ order 


NNLL 


3-loop 


2-loop 


1-loop 


LO 


3'"^ order 


N^LL 


4- loop 


3-loop 


2-loop 


NLO 


4**" order 


N^LL 


4-loop 


3-loop 


3-loop 


NNLO 



The first three orders correspond to traditional counting in renormalization-group improved 
perturbation theory, while 4*'^order simply uses all the available information. 

Next, we consider, the separate variations. The bands in the first four panels of Figure [7] 
show the effect of the scale uncertainties. The bottom two panels of Figure [7] show the effect 
of the more natural correlated and anti-correlated scale variations introduced in [11]. The 
correlated variation is defined to hold fJ^j/fJ^s fixed. So we vary 



/ij -> c^/TQ, fis -^ ctQ , - < c < 2 . 



(41) 



This probes the upper and lower limits on fij and fig, but avoids the unphysical region where 
fig < fJ'j or fih < jJ'j- The orthogonal anti-correlated variation is defined to hold fj,'j/{Qfis) 
fixed. It is 

/i| -)■ aQ'^T fis -^ aQr, —^<a<\Pl. (42) 

v2 

This is independent from the correlated mode but again avoids unphysical scale choices. 

Overall, we find good convergence order-by-order in perturbation theory. However, some of 
the higher-order scale variations are outside of the range of the lower orders. This was not the 
case for thrust, where the central value of the prediction was much more stable. Nevertheless, 
for both thrust and heavy jet mass, the complete perturbative uncertainty, defined as the 
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Figure 7: Perturbative uncertainty at Q = 91.2 GeV. Each of the scales is varied separately 
by a factor of two around the default value. We show ^p^ and, for reference, aleph data at 
LEP 1 scaled by the central value of each bin. All plots have as{mz) = 0.1214. 



envelope of the various variations (that is, the maximum and minimum over them) does have 
the higher-order bands contained within the lower order bands, as can be seen in Figure El 

Next, we fit the theoretical prediction to the aleph data from 91.2 to 206 GeV [23] and 
extract a^. The fit is done by minimizing the x^^ using experimental statistical uncertainties, 
for the theory prediction with default scale choices. The statistical error on a^ is determined 
by variations around this minimum. The perturbative uncertainty is extracted with the 
uncertainty band method [31], exactly as in [TT] for thrust. The envelope over the hard, 
matching, correlated and anti-correlated scale variations are included in this extraction. We 
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Q 


91.2 


133 


161 


172 


183 


189 


200 


206 


AVG 


fit range 


0.08 

0.18 


0.06 
0.25 


0.06 

0.25 


0.06 
0.25 


0.06 
0.25 


0.04 
0.20 


0.04 
0.20 


0.04 
0.20 


- 


xVd.o.f. 
stat. err. 
syst. err. 
hadr. err. 
pert. err. 
soft. err. 
tot. err. 


67/9 
0.0002 
0.0011 
0.0044 

+0.0006 
-0.0011 

0.0005 
0.0047 


2.3/4 
0.0055 
0.0011 
0.0028 

+0.0006 
-0.0011 

0.0005 
0.0064 


0.66/4 
0.0108 
0.0011 
0.0022 

+0.0009 
-0.0014 

0.0006 
0.0112 


1.8/4 
0.0144 
0.0011 
0.0021 

+0.0003 
-0.0005 

0.0002 
0.0147 


5.2/4 
0.0065 
0.0012 
0.0019 

+0.0007 
-0.0011 

0.0005 
0.0070 


1.1/4 
0.0032 
0.0013 
0.0018 

+0.0006 
-0.0009 

0.0004 
0.0040 


8.8/4 
0.0034 
0.0014 
0.0017 

+0.0006 
-0.0008 

0.0004 
0.0041 


3.8/4 
0.0034 
0.0011 
0.0016 

+0.0005 
-0.0007 

0.0004 
0.0040 


0.0015 
0.0013 
0.0022 
0.0009 
0.0004 
0.0031 


as{mz) 


0.1214 


0.1235 


0.1328 


0.1077 


0.1267 


0.1234 


0.1218 


0.1189 


0.1220 


PYTHIA 
ARIADNE 


0.1365 
0.1238 


0.1239 
0.1262 


0.1333 

0.1355 


0.1073 
0.1093 


0.1266 
0.1288 


0.1214 
0.1239 


0.1202 
0.1731 


0.1168 
0.1687 


0.1230 
0.1250 



Table 1: Best fit to aleph data. The row labelled "pert err." is derived from scale uncertain- 
ties and the row labelled "soft err." from the uncertainty on cf and cf^ in Eq.( l33l) . The rows 
labeled pythia and ARIADNE give the value of a^ after correcting for hadronization and quark 
masses using pythia or ARIADNE. The ARIADNE corrected prediction for the two highest two 
energies produce very poor fits, and are excluded from the average 



also include an additional soft uncertainty associated with the errors in the extraction of cf and 
cf^- These are computed by fitting as within the errors on cf and cf^ in Eq. (!33|l . and taking 
the difference with the central value as the uncertainty. The soft and perturbative uncertainties 
are assumed uncorrelated. The systematic uncertainties are taken from [18]. To use these 
uncertainties, we are forced to keep our fit ranges the same as in ^E\. The hadronization 
uncertainties are also taken from [TS], which are based on Monte Carlo simulations. Note 
that, as in [TT], we use the uncertainties from [TH] but do not correct for hadronization. 
Hadronization will be discussed in detail in Section |5l Finally, the values for each energy 
are combined with a weight inversely proportional to the square of that energy's total error. 
The statistical uncertainties are assumed uncorrelated, and combined in quadrature, while for 
the other uncertainties a linear weighted average is performed. The results are tabulated in 
Table [H 

We show in Figure [8] the convergence of the best fit values as a function of energy. There 
is very good consistency among the different energies and the convergence order-by-order in 
perturbation theory is good as well. The fit values for different orders are given in Table [2l 

The final fit for heavy jet mass gives 

as{mz) = 0.1220 ± 0.0014 (stat) ± 0.0013 (syst) ± 0.0022 (had) ± 0.0009 (pert) ± 0.0004 (soft) 
= 0.1220 ± 0.0031 (Heavy Jet Mass) . (43) 
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Figure 8: Best fit values for as{mz)- From right to left the lines are the total error bars at 
each energy for first order, second order, third order and fourth order, as defined in the text. 
The bands are weighted averages with errors combined from all energies. 





LEP 1 +LEP 2 


LEP 1 (91.2 GeVonly) 


order 


as 


total err 


pert, err 


as 


tot. err 


pert, err 


P*order 
2^'^OTdeT 
3'''^order 
4*^order 


0.1111 
0.1156 
0.1189 
0.1220 


0.0104 
0.0064 
0.0038 
0.0031 


0.0100 
0.0057 
0.0025 
0.0009 


1099 
0.1132 
0.1168 
0.1214 


0.0100 
0.0072 
0.0052 
0.0047 


0.0110 
0.0055 
0.0026 
0.0011 



Table 2: Best fit values and uncertainties at different orders. 

This can be compared to the result for thrust, using exactly the same technique, and the same 
energy ALEPH data (Table 2 of [H]). Updating this result to include the more recent NNLO 
distributions [3111], using the same cf values, Eq.( l33|) . with associated "soft" uncertainty, and 
restricting to only the ALEPH data, we find 

asimz) = 0.1175 ± 0.0009 (stat) ± 0.0011 (syst) ± 0.0014 (had) ± 0.0016 (pert) ± 0.0006 (soft) 
= 0.1175 ±0.0026 (Thrust). (44) 

Combining these results, assuming 100% correlation between heavy jet mass and thrust, gives 

as{mz) = 0.1193 ± 0.0011 (stat) ± 0.0012 (syst) ± 0.0017 (had) ± 0.0013 (pert) ± 0.0005 (soft) 
= 0.1193 ± 0.0027 (Combined) . (45) 

This value is consistent with the recent world average of asirriz) = 0.1184 ± 0.0007 [16j. 
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5 Non-perturbative effects and quark mass corrections 



The Us fit from the previous section used the theory prediction at the parton level with five 
fiavors of massless quarks, neglecting hadronization and quark masses. Hadronization induces 
a power correction on the heavy jet mass distribution. Its effect is suppressed by a small scale, 
such as Aqcd/Q oi' ^qcd/ l^s relative to the perturbatively calculable part. The 6-quark mass 
corrections are suppressed by m^jQ. These effects are therefore both parametrically smaller 
than the large logarithmic corrections which we resum. Nevertheless, they are quantitatively 
important, and our final uncertainty on a^ is dominated by the way these power corrections are 
modeled. The dominant part of the 6-quark mass corrections is calculable, and is expected to 
shift as at around the 1% level, as observed in [II1[T1]. The inclusion of 6-quark corrections will 
be an important addition for future work. However, since they scale like l/Q, the dominant 
effect of these mass corrections can be absorbed into the same power correction model as 
hadronization effects, which also scale as l/Q. In this section, we explore the Monte Carlo 
treatment of power corrections, and an alternative theoretical model. 

Monte Carlo simulations can include quark masses explicitly. They also attempt to model 
hadronization, for example with a string fragmentation model in pythia. This produces an 
event with stable particles which can be run through a detector simulation. Such simulations 
are an essential part of every experimental study, and must play some role even for inclusive 
event shape analysis. For example, the event shape is often measured using only the charged 
particles, whose momenta are more precisely known, and then corrected to all particles with 
help of the simulation. Monte Carlo hadronization models have a number of free parameters 
and can usually be tuned to any particular data set so that the simulation reproduces the 
data quite well. However, no single tuning reproduces all the data, and therefore different 
tunings are often used for different analyses. A more troubling fact is that, as demonstrated 
in [TT], the tunings often correct for features having nothing to do with hadronization, such as 
subleading log resummation. Such tunings are guaranteed not to scale well with energy. This 
may be a serious problem for high energy colliders which simultaneously probe many energy 
scales, such as the Large Hadron Collider at CERN. 

The hadronization uncertainty we used for the a^ determination were taken from [TB] , but 
we have also studied hadronization and mass effects in the Monte Carlos on our own. The 
last two rows of Table [T] show the best fit vales for a^ after the theory is corrected bin-by- 
bin for both hadronization and bottom and charm quark mass corrections using the Monte 
Carlo event generators pythia v6.409, with default parameters [23], and ARIADNE v4.12 
with the ALEPH tune [32]. Recall that ARIADNE actually feeds through pythia to handle 
hadronization, so the difference is entirely due to the way the parton shower is implemented. 
With thrust, the same exercise was performed, and the corrections with ariadne were found 
to be very small, which helped justify not correcting for hadronization and quark masses at all 
in the published as fits. For thrust corrected with pythia, there was a systematic downward 
shift in as- For heavy jet mass, the corrections with ariadne are large. In fact, for the 
high energy data, ridiculous values such as as = 0.1731 result. The pythia corrections are, 
for heavy jet mass, smaller than they are for thrust. In fact, we find a bigger discrepancy 
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Figure 9: Hadronization and mass corrections with PYTHIA. The theoretical prediction using 
PYTHIA at the hadron level with massive quarks and the parton level with massless quarks is 
compared to data and to the 4*^order theoretical prediction using SCET. For thrust, pythia 
agrees remarkably well with the data, while for heavy jet mass, there is a substantial discrep- 
ancy especially in the fit region, which is zoomed in on in the bottom panels. 



between the thrust and heavy jet mass Ug fits after correcting with either Monte Carlo than 
without. Thus, although we cannot justify correcting the theory curve with either Monte 
Carlo, we confirm that the hadronization uncertainties listed in Table [H which were taken 
from [18], span reasonable Monte-Carlo simulated variations due to hadronization and quark 
mass effects. 

To understand why the power corrections come out so differently for thrust and heavy jet 
mass, we compare PYTHIA at the parton and hadron levels to the 4*'^order SCET prediction 
(N^LL + NNLO), and to the aleph data at 91.2 GeV in Figured From the top two panels, 
we see that in the peak region, in both cases the parton-level theory prediction comes out 
somewhere between the parton and hadron level Monte Carlo. However, in the bottom two 
panels, which zoom in near the fit region, the difference between the two event shapes is much 
more dramatic. For heavy jet mass, the SCET curve is above the data, while partonic pythia 
is below it and hadronic pythia is even farther below. In contrast, for thrust, all of the curves 
are much closer and the power corrections, as modeled by pythia are a much smaller effect. 
It is clear that pythia has trouble handling both event shapes simultaneously. 

An alternative to using Monte Carlo simulations to simulate hadronization is to model 
the power corrections directly with effective field theory. As discussed in [22], hadronization 
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effects can be absorbed into tlie soft function by convolution of the perturbatively calculable 
part with a non-perturbative shape function 

Sin\\{kL, kn, /i) = / dk'^dk'j^SpartikL - ^l, kn - /c^,/i )5'mod(fcL, k'j^) , (46) 

where S'part(^L, k^, /i) is what we have previously just been calling S{kL, kn, jj) and S^odiki, kjij 
is a non-perturbative model function. Generally, S^odiki, kjij is expected to have support only 
for ki, kji < Aqcd- As observed in [33], there is an ambiguity in the factorization of the soft 
function into these two pieces, which leads to a difficulty in assigning physical significance 
to SmodikLykn) and poor convergence in perturbation theory. This ambiguity is associated 
with the existence of a renormalon, which can be removed within SCET [331 122] • Indeed, if 
data closer to the peak region were included in the fit, or if convergence of the model function 
parameters were an issue, removing the renormalon could have an important effect. Since 
we are not immediately interested in these issues, for simplicity, we will simply ignore the 
renormalon. 

The simplest model function is just composed of delta functions 

Smod{kL, kn) = 6{kL - -Anp)S{kR - -Anp) • (47) 

The one parameter, A^p, can be thought of as representing the mass gap of QCD due to 
hadronization and therefore should be of order Aqcd- This model function allows us to fit the 
leading power correction. Any other one-parameter family of shape functions can be written 
in this form up to higher power corrections, which should have a subleading effect on the 
distributions. For example, the smallest scale probed in our fits is the soft scale at the lower 
end of the fit region at 91.2 GeV, ^ = ^, = pQ > (0.08) (91.2 GeV) ~ 7 GeV. With Aqcd ~ 300 
MeV, this can be a 4% effect. Higher power corrections, of order (Aqcd//^)^ should have less 
than a 0.2% effect in our fit range. 

Once this shape function is convoluted with the perturbative distribution, it has the effect 
of simply shifting the distributions 

da, , dcr, . . .,„, 

-(r)^-(r-ANp) (48) 

dcr , , dcr , 1. , , , , 

^(P)-^(P-2A,P). (49) 

The factor of | is easy to understand. The shift causes each hemisphere mass to increase by 
^Anp. Since thrust sums both hemisphere masses, while heavy jet mass measures only one, 
heavy jet mass feels only half of the increase. 

This model was studied for thrust in [TT], where it was found that a larger Anp can be 
compensated for by smaller as leading to a flat direction in the two parameter flts. We 
reproduce this result in Figure [IDl This flgure shows the 2(T and 5cr confldence regions in a 
combined flt to all of the aleph data for thrust from LEP. On the same plot, using the same 
model function, we show the contours for heavy jet mass. First of all, we observe that the 
flat direction exists in both of the data sets. We might have hoped that having two event 
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Figure 10: Contours of 2a and 5a confidence in the simultaneous fit of a^ and a non- 
perturbative shift parameter Anp to the thrust and heavy jet mass aleph data from 91.2 
to 206 GeV. The combined fit is also shown. 



Event Shape 


as{mz) 


Anp (GeV) 


xVd.o.f. 


Thrust 

Heavy Jet Mass 

Combined 


0.1101 
0.1017 
0.1236 


0.821 

3.17 

-0.621 


66.9/47 
60.4/43 
453/92 



Table 3: Best fit values including leading power correction. The x^ is calculated using both 
statistical and experimental systematic uncertainties. 



shapes would remove the ambiguity, but this does not happen. Second, we see that while the 
perturbative fit has a^ lower for thrust than for heavy jet mass, with the power corrections, 
the value of a^ is higher for thrust, as found in previous studies pTf \TE\ . However, when we 
perform a simultaneous fit to all of the thrust and heavy jet mass degrees of freedom, we get a 
value for as that is larger than each one separately. The best fit for thrust, heavy jet mass, and 
the combined fit are shown in Table [31 The fact that the thrust and heavy jet mass contours 
do not overlap indicate that a better handling of non-perturbative effects is required. 

We conclude that neither correcting the theory curves with a Monte Carlo simulation nor 
using a minimal shape function approach for the leading power correction is satisfactory. The 
shape function approach is improvable, while the Monte Carlo approach is limited by the 
perturbative accuracy of the parton shower, which will be limited to leading-log resummation 
in at least the near future (although SCET may eventually help go beyond LL [SU |35]). To 
improve the shape function fit, a number of additional ingredients should be included. First 
of all, the renormalon ambiguity in separating the perturbative and non-perturbative parts of 
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the soft function should be removed. This is not hkely to have much effect in the fit region 
we use, but will allow us to use data closer to the peak. Having more data involved will more 
highly constrain the fit and could remove the flat direction. To do this, we would need the bin- 
by-bin correlations among the experimental systematic uncertainties, which are not publicly 
available. In addition, there are perturbatively calculable effects we have not included, such 
as electroweak and mi, corrections, as in [T3], which may have up to a 1% effect. It would be 
very interesting to see if the thrust and heavy jet mass distributions can be reconciled once a 
thorough effective field theory analysis, including non-perturbative effects, is performed. 

6 Conclusions 

In this paper. We have studied the heavy jet mass distribution using Soft-CoUinear Effective 
Theory including N^LL resummation and matching to the NNLO fixed order distribution. Up 
to this point, this kind of accuracy has only been available for the thrust distribution. Having 
an additional event shape helps control for systematic uncertainties, making the fit for a^ 
more trustworthy. It also gives us insight into power corrections and multi-scale soft functions 
which will be important for the LHC 

The heavy jet mass fit gives asirnz) = 0.1220 ±0.0031. This value is larger than what had 
been found for thrust at the same accuracy, asimz) = 0.1175 ± 0.0026. The uncertainty on 
heavy jet mass is larger partially due to a larger hadronization uncertainty. In our study, no 
corrections were made for hadronization. We explored the traditional method of hadroniza- 
tion, using Monte Carlo event generator, such as pythia and concluded that such an approach 
is problematic for theoretical calculations of this accuracy. Since the Monte Carlo has been 
already tuned to the data we are trying to match, the tuning has partially compensated for 
resummation of subleading logarithms. Comparing pythia's output in the fit region, the 
hadronized distribution is actually farther away from the data than the parton-level distribu- 
tion. 

Our as values from thrust and heavy jet mass contrast with the results of [T7], which at 
NLL+NNLO accuracy derived as{mz) = 0.1266 from thrust and as{mz) = 0.1211 from heavy 
jet mass. A comparison of various fits to thrust and heavy jet mass using the same aleph 
data is shown in Table HI The authors of [17] have observed that event shapes tend to belong 
to one of two classes. The first class, including thrust, tends to produce higher values of a^ 
than the second class, which includes heavy jet mass. These authors attributed the difference 
to better perturbative stability in the second class. We find, if anything, better perturbative 
stability for thrust. Instead, the reason for the systematic separation of as values in this study, 
and also in the NNLO study of [18], may have more to do with their use of a Monte Carlo 
simulation to correct for hadronization. A similar conclusion was reached in [3S] which studied 
event shape moments. The values of as for the two classes must eventually be reconcilable, 
but there may be a physical reason why the power corrections for one class are larger than for 
the other. This is worth understanding more thoroughly, and may have implications for the 
LHC. 

The alternative to using a Monte Carlo simulation for hadronization is to add a shape 
function contribution within the effective field theory. Our simple shape function study shows 
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Order 


N^LL+NNLO 


N^LL+NLO 


NNLO 


NNLO n] 


NLL+NNLO [n] 


hadronization 


NO 


NO 


NO 


YES 


YES 


Thrust 
Heavy Jet Mass 


0.1175 
0.1220 


0.1173 
0.1189 


0.1262 
0.1265 


0.1275 
0.1248 


0.1266 
0.1211 



Table 4: Best fit values for as{mz) at various orders in perturbation theory. The first three 
columns are our results, the last two which include a Monte Carlo based hadronization cor- 
rection are from |17] . 



that the leading power tends to shift a^ from both heavy jet mass and from thrust to lower 
values, with the heavy jet mass shift of larger magnitude. This can help explain why the thrust 
as comes out lower than the heavy jet mass a^ in our study, and not in [T71 [18]. However, 
we also found that the best fit over all the aleph data from 91.2 to 206 GeV for thrust was 
incompatible with the best fit from heavy jet mass, and that the fiat direction between as and 
the non-perturbative parameter Anp persists in both distributions. 

To get the values of as extracted from thrust and heavy jet mass to agree may require 
including additional ingredients, which can be done within the effective field theory framework. 
For example, there is a calculable mh correction which tends to bring as up at least for 
thrust [11] . Including every possible correction must produce the same value of as from thrust 
and heavy jet mass, and it will be interesting to see precisely how this happens. Also, more 
data should be included. Using data for values of heavy jet mass and thrust closer to the peak 
will lead to a more constrained shape function fit, although it may require going beyond the 
leading power. In addition, using data from other lep experiments and other experiments at 
lower center-of-mass energy can further test and constrain the event shapes. 

However, it is not clear if all of the differences between thrust and heavy jet mass can be 
accounted for entirely within SCET. For example, there is the possibility that the difference 
between thrust and heavy jet mass has more to do with the way hadron masses are handled 
experimentally than from higher order power corrections. In [37], substantial differences in the 
form of power corrections among the ii^-scheme, p-scheme and decay-scheme were found. It 
may turn out that an ultra-precise as fit can only be made if the identity of all the hadrons is 
known, which may be possible for future measurements but is not available for existing data. 
More likely, the thrust and heavy jet mass distributions can be made to agree within SCET, 
but the uncertainty on a^ will ultimately be limited by a hadron-mass-scheme dependent 
uncertainty. In any case, once the ingredients discussed for thrust in [H] are applied to heavy 
jet mass, we will be able to extract a more precise lesson about the importance of power 
corrections. In addition to reducing the uncertainty from as and teaching us about power 
corrections, combining the insights from thrust and heavy jet mass will more generally pave 
the way for deeper understanding of relevant jet-based observables at the LHC. 
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Figure 11: A comparison of the full fixed-order calculations and expanded SCET at NLO. 
Update of Figure [1] with cutoff y = 10~^^ in B{p) from event 2. 

Note Added 

After this manuscript appeared, it was suggested that the precision on cf and cf could be 
improved by lowering the infrared cutoff used by EVENT 2. The cutoff y is implemented by 
throwing out events if two partons have {pi + ^2)^ < yQ"^- The default cutoff is 10~®, and 
the authors of event 2 caution about numerical instabilities if the cutoff is taken too low. 
We find that for cutoffs below 10~^^, there are insurmountable numerical problems, however 
y = 10~^^ seems to be convergent. We therefore ran 135 billion events with y = 10~^^ and 
2500 bins (Ap = 0.0002)) - the main text uses 10 billion events with y = 10~^ and 1000 bins 
(Ap = 0.0005). The difference between this new numerical data and the SCET prediction for 
the singular terms is shown in Figure [TTl which is to be compared to Figure [H One can see 
that the curves for all color structures now converge to zero, as expected. 

Next, the constant cf was extracted from these curves. The value cf for various lower 
bounds pmin are shown in Figure [121 Again, improved numerical stability is clear. Fitting 
a sixth order polynomial to the 59 points between p^^^ = 0.0004 and pmin = 0.012 and 
extrapolating to pmin = leads to 



-ap 



(49.1)C^ + (-33.2)Ci.C^ + (50.2)C^Ti.n/ 



(50) 



The value of the Cp coefficient is now consistent with the prediction of ^ = 48.7 from non- 
Abelian exponentiation. The fit is somewhat sensitive to the lower value of pmin used in the 
regression, but not very sensitive to the upper value. Fitting a fourth order polynomial to 
the 38 points between p^in = 0.0006 and p^m = 0.008 gives cf^ = (49.8)^2 + (-33.3)C^C^ + 
(50.3)Ci?Tpn/. Since the CfCa and CpTprif terms are practically unchanged, and the Cp 
term is fixed by non-Abelian exponentiation, it is reasonable to assume that the remaining 
uncertainty on these numbers will have a negligible effect on the a^ fits. Performing the same 
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Figure 12: Extraction of the two-loop constants in the soft function. Update of Figure |2] with 
cutoff y = 10"^^ in B{p) from EVENT 2. 



analysis for thrust, leads to 



cf = (49.1)^2 + (-57.8)C^C^ + (43.4)C^Tirn/ 



(51) 



Combining these, assuming the Hoang-Kluth Ansatz for the soft function, Eq. ( 130|) . gives 

cIl = (0)C| + {-7.5)CpCa + {-2.1)CFTFnf (52) 



Thus, we take 



71 



r^ 

'-2L 



^Ci 



io)cl 



-57.8)CfCa + {'i3A)CFTpnf 
-7.5)CfCa + (-2.1)CpTpn/ 



(53) 
(54) 



With these more accurate numbers and a more accurate numerical calculation of the NNLO 
distribution, we can now repeat our comparison of the singular terms to the exact distribution. 
Using an infrared cutoff of 10~^ for the C functions, the agreement with the singular terms is 
improved. This can be seen in Figure [T^ which is an update of Figure [31 Taking the difference 
between the curves gives Figure [HI One sees that the 1/N^ color structure, corresponding to 
Cp, has improved convergence towards zero. If these curves were known with perfect accuracy, 
they could be used to test the Ansatz in Eq. (j30|) . The most poorly convergent color structures, 
1/N'^ and n'j are not sensitive to this Ansatz, and the others are consistent with convergence 
to zero within the statistical uncertainty on the numerical NNLO calculation. 

Finally, we reconsider the «« fits in light of these more precise soft function coefficients and 
NLO matching functions. Refitting the thrust distribution to the aleph data changes asirriz) 
from 0.1175 to 0.1176 and refitting the heavy jet mass distribution raises asirriz) from 0.1220 
to 0.1224. These shifts are within the quoted soft uncertainties. 
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Figure 13: Comparison of tlie full NNLO heavy jet mass distribution and the singular terms. 
Update of Fig [3] with cutoff y^ = 10^^ in C{p). 
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A Soft function 

To evaluate the heavy jet mass distribution with NNLO precision, we need the soft function 
at the scale fig evaluated to order a^. The Laplace transformed soft function can be written 
as 

s{Li,L2) = s^(Li,/i)s)j(L2,/i)s/(Li -L2), (55) 

where Li = ln{fiULe'^'^) and L2 = ln{ fiUfie'"^), with ul and u^ the Laplace conjugate variables 
to the soft momenta ki and fc/j. The /i-dependence is determined by the function we call 
's^{L,fi). This is equivalent to the function Us{x, fi, {ixe^'^)~^) in [22], and we have already 
given its expansion to order a^ in Eq. f lT^ . The function 'sf{L) is //-independent with a^ 
evaluated at the scale {i^li^r)~^^'^- It is more useful to be able to use a^ evaluated at the scale 
/i, which we can do with the replacement 



as^ fas{iJ,)\ , fasifJ^y 



AnJ ■ V 4- ; ' V 4- ; Mo(^i + ^2)] + (^)'[/3o^(L, + L2)^-/3,(L, + L2)]. 

(56) 
Then truncating the soft function to finite order will induce some residual /i- dependence. 
For the finite part 'sf{L), we use for numerical studies the form in Eq. (130|) : 

ML) = 1 + (^) <^f + i^f ^f + 4lL'] +■■■. (57) 



with 
and 



-Cpn^ (5^ 



TT^ 



—C'p + {-57.8)CfCa + {A3A)CpTFnf (59) 



4l = (O)C^ + i-7.5)CFCA + i-2.1)CFTFnf (60) 

which have been extracted using SCET from the exact NLO thrust and heavy jet mass distri- 
butions. (Note: these numbers are updated to Eq. (!53|) . The main text and fits use the earlier 
values in Eq. ( 15^ .) 

More generally, for NLO-matching, all that is relevant is a single projection of the order 

a^ soft function 

1 



c 



s 

2p 



TT 



Sf2mdL = 4--cl^. (61) 
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For NNLO matching, another projection is necessary, of the form in Eq. f l27p . For the form 
in Eq. (E7]), this projection is 



-2C 



TX 



Sf2{iL)\n 



2cos( 



,L^ 



dL = 4C34i 



These are then expanded as 



(62) 



(63) 



C2 = CpC2Cp + CfCaC2Ca + CFnFTFC2nf ■ 

We will use these expressions for the singular heavy jet mass expansion and the Gij coefficients 
below. We also use anomalous dimensions and /3-function coefficients which can be found 
in [IT]. 



B Expanded soft function 

Putting the pieces together, the soft function expanded to order a^ with a^ = ^^(/i) is 



s{L,, L2, /i) = 1 + (^^ j [ - (Li^ + L2')ro + (Li + L2h^ + cf 

+ {^y lliLi' + L2')Tl + L,L2{L,To - 7o^)(I^2ro - 7o^) + (^1' + ^2') Q/^oFo - ToTo^ 

+ {L^ + L2') (-T, - /3o7o^ + 1^0 f - cfro) + {Li + L2)(7f - cf/3o + cf7o^) 
+ cf + (Li - L2)'cfi 



+ -7o^r2(Li% + L2%) 



3r 2\ 



/3or^ + 75'r^ (Li^V + VLi 



2 r 2 

2 



+ ( - ^(7o^)'ro + li'o^oTo - l^'oTo + rori)(L/ + L2') 

+ ( - (7o^)2ro + ^7o^/3oro) (^^1=^1^2 + ^2'^i) + ( - (7o^)'ro + 2j^f3oTo + 2ToT,y,^L 

+ (^(7o^)' - ilo)'(3o + ^7o^/3o' - 7f ro + ^/^iFq - 7o^ri + ^/3oFi) (L^^ + 1.2=^) 

+ ( + ^(7o^)' - (7o^)'/3o - 7f ro - 7o^ri) {L^L2 + L2^L,) 

+ (7o^7f - 27f /3o - 7o^/3i - ^2) {U^ + 1^2') + 2^hiLiL2 + 72''(^i + ^2) 

+ cf |^F2(Li^ + 2Li2L2' + L2^) + ( - 7o'ro + ^-P^T^Yu^ + ^2^) 
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+ ((7o^)^ - 2^1 ^.y^L^ +(-/?! + 7f) (^1 + ^2)} 



+ [4+4l{Li-L^ 



-2/3o + 7o^)(^i + ^2)-ro(Li2 + L2' 



(64) 



C Singular terms in the heavy jet mass distribution 

The singular part of the heavy jet mass distribution is calculable in SCET. Writing it in the 
form 



«s^3 



D{p) = 5{p)Ds + (^) [DM]^ + (^)'[^b(p)]+ + [^] [Dc{p)V + ■■■ 



the result is 

— (S 



Cf\-2 + 



27r^ 



+ 



a.s 



An 
+ CpTpUf 



Cj(4 + ^-48C3)+C.C.(^+«^'^ ^^'*'^«««^ 



vr 

To " 

28 147r2 
81 3~ 



6 



90 



9 



9 l^^^'P 



and 
DAip) 

Db{p) 



C, 



cl 



-41np-3 } , 
81n3p + 181nV+(13 



— )lnp + --7r2-4C3 



22 

+ CpTpUf - 4 InV H In p + 5 

9 



+ CfCa 



1 



11 InV + ( + ) Inp + 6C3 

'^ ^ 18 3 ^ '^ 4 ^'^. 



Dc{p) =-\Cl -81n^p-301nV + ln^p( -44 + 87r2) +lnV(8C3 + 127r^-27 



+ lnpf -c; 



s 

2pCf 



48C3 



4l7r4 137r2 17 



47r2 
+ ^-C3 + 14C3 + I2C5 



90 3 2 

3tt* Svr^ 47 3 



40 4 



g 4'^2pCF 2'^2CCf 



(65) 

(66) 
(67) 
(68) 



(69) 



+C'pnfTp 



401n% 561n^p , , 

H ^ + In^ p 

3 9 '^ 



43 



167r2 



lnp(^ 



/232C3 587r2 1495 



9 



9 



^2pnf 



29 



254C3 Tvr^ 71 vr 



9 15 

CpnJT} 



+ 



1511 

los" 






1121n^p 681nV , /140 IGtt^ 



176C3 647r2 



27 



3598 2 5 - 
^43" ^ 3^2""^. 



-C'fC'a 



8471nV , 2 /3197 llvr^ 
'^V 36 3 



27 



, , , llvr^ 857r2 11323\ 

Inp 22C3 + 

^^ ^^ 45 9 324 / 



, 36IC3 54l7r^ 



7397r2 



77099 11 



+CpCA 



27 ■ 540 
110 InV 



+ ln^p(-y 



486 



57r' 



Q ^2pCA 



+ ln> -36C3 + 



357r2 467 



+ Inp 



1682C3 1337r^ 4037r^ 
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-3OC5 



1943C3 27r2C3 777r4 75771^ 49 11 



+CACFnfTF 
, 6O8C3 



+ 



616 



40 



72 



29663 
324 

_ 1 5 

2j Q ^2pCp ^^2pCA 'f-^C^^A 



-2pCA 



fAii^ 



lnV + ln>(^— 



512\ 



Inp 8C; 
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D G^o expansion 



Occasionally it is helpful to write an event shape distribution as 

CO i+1 



R{x) 



V m=l / \i=l 7=1 / n=0 



[X 



(71) 



The Gij and Cm are calculable in SCET for exponentiation up to N'^LL accuracy. 
The results are 



---K-I^t)' 



^ V 8 40 2 ^^4 ^''^^ 



+ C^Cf( 



648 



360 
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G12 — —2Cf 
Gil = 3Gi7' , 



CpUfTp ( 



905 58 



V162 9 



■Cs 



77r2 1 . 
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8977 737r^ 8571^ 481 
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^23 = Cp 
Ci'22 = (^F 
G21 = Cf 
G34 = Cp 
G33 = Cf 

CpUfTp [ 2 



^4 111 



(73) 



-C^^ + n/Tpy + C^ 



169 7r2 

y ^y 



C^(4C3 + ^) - ^n^Tp + Ca(^ - 6C3 



108 ^ -^ ^ 27 



ATTT^ ^ C'a%T^1^ - ^fTpi^^ 
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27 
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f ^27 



^ ^ llTT^ ^2l6 ■ 
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77 77r^ llTT^ 18i 
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